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Our formalism described recently in [13] is applied to the study of particle 
creation in spatially uniform electric fields, concentrating on the cases of a 
time-invariant electric field and a so-called 'adiabatic' electric field. Several 
problems are resolved by incorporating the 'Bogoliubov coefficient' approach 
and the 'tunnelling' approaches into a single consistent, gauge invariant formu- 
lation. The value of a time-dependent particle interpretation is demonstrated 
by presenting a coherent account of the time-development of the particle cre- 
ation process, in which the particles are created with small momentum (in the 
frame of the electric field) and are then accelerated by the electric field to make 
up the 'bulge' of created particles predicted by asymptotic calculations [2]. An 
initial state comprising one particle is also considered, and its evolution is de- 
scribed as being the sum of two contributions: the 'sea of current' produced 
by the evolved vacuum, and the extra current arising from the initial particle 
state. 
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1. INTRODUCTION 

It has long been known that a strong classical electromagnetic back- 
ground can cause vacuum instability and generate particle/antiparticle 
pairs. This phenomenon was first studied in 1951 by Schwinger [1], who 
predicted particle creation in a constant electric field. Schwinger's pre- 
dictions were based on the imaginary component of the effective action, 
but the same result has since been derived from a variety of alternative 
approaches: Bogoliubov coefficients within a Canonical formalism [2], tun- 
nelling amplitudes [3, 4, 5], a wave- functional method [6, 7], and particle 
detectors [8, 9]. (In [10, 11] back-reaction is also considered.) A brief 
review is presented by Sriramkumar et al. [8], who point out that the 
predictions of these approaches do not always agree. The same authors 
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point out in [12] that approaches based on Bogohubov coefficients or on 
tunnelUng ampUtudes are Hable to suffer from a gauge-dependent particle 
definition (i.e., choice of 'in' and 'out' modes) which can in turn lead to 
gauge-dependent predictions of the number of particles created. 

We have recently introduced an alternative approach to the study of 
particle creation in electromagnetic backgrounds [13] which uses methods 
analogous to those of conventional multiparticle quantum mechanics. Our 
approach is to emphasise the states of the system, described in terms of 
Slater determinants of Dirac states. The vacuum state 'at time t' is de- 
scribed in terms of the Slater determinant of all 'negative energy' states, 
and provides a concrete realisation of the Dirac Sea. This approach has 
provided straightforward derivations of the general S-Matrix element and 
expectation value of the theory, using only methods familiar from mul- 
tiparticle quantum mechanics. We argued in [13] that this approach is 
often quicker and clearer than conventional methods. More importantly, 
the 'Bogoliubov coefficient' and 'tunnelling amplitude' methods were com- 
bined consistently into a single gauge invariant formulation, which could 
be linked to the motion of an observer (detector) . 

As with the Canonical, Tunnelling, and Wave-functional approaches, our 
analysis depends explicitly on a choice of 'in' and 'out' particle interpreta- 
tion, and on the corresponding categorisation of in and out modes. Var- 
ious categorisation schemes have been proposed, based on asymptotic or 
adiabatic properties of solutions, or on the diagonalisation of a suitable 
Hamiltonian. Unfortunately, most schemes depend on a choice of coordi- 
nates (as in the Unruh effect [20]) or of gauge [12], and the relation of these 
to the behaviour of a particle detector is often ambiguous [8]. Also, those 
schemes based on asymptotic or adiabatic approximations generate accu- 
rate predictions only at asymptotically late times or in sufficiently weak 
electromagnetic backgrounds. Particle detectors provide a more opera- 
tional particle concept, but their predictions are not always proportional 
to the number of particles present [8, 9] even when the detector is inertial. 

The categorisation scheme suggested in [13] provides a consistent parti- 
cle interpretation at all times without requiring any asymptotic conditions 
on the in and out states. It consistently combines the conventional 'Bo- 
goliubov coefficient' and 'tunnelling amplitude' methods, resolving gauge 
inconsistencies that trouble each of these [12]. Also, by utilising the con- 
cept of radar time, this definition naturally incorporates the motion of the 
observer (detector), providing a definition which depends only on the ob- 
server's motion and on the background, and not on the choice of coordinates 
or gauge. 

The present paper presents a simple application of this formulation. We 
shall examine the problem of pair creation in time-varying spatially uni- 
form electric fields, as viewed by an observer stationary with respect to 
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this electric field. We present definite results for a constant electric field 
E{t) = (0,0, E) and for the so-called 'adiabatic' [2] electric field E{t) = 
(0,0,-Escch^ (mt//?)) . Such backgrounds have been studied by various au- 
thors [1, 7, 5, 4], who derived pair-creation rates as time-averages over an 
infinite period of time. The problem of pair creation after a finite evolution 
time has only recently been considered only recently [6, 2]. The electric 
field was taken as constant within a certain time interval (say 0<t <T), 
but was zero before and after. The 'in' and 'out' states were defined in the 
'free regions' before the field was turned on (t < 0) , and after it was turned 
off {t>T), where the solutions were simple plane wave states. Dependence 
of the results on T was then studied. 

We use the particle definition proposed in [13] to consider the backgroimd 
E= {0,0, E) for all time, and derive the magnitude of pair creation in a 
state which has evolved from the vacuum for a finite length of time. That 
is, we answer the question, "Suppose the system is in the vaciium state at 
some time to. What will be the properties of the system at time to + T"?" 
We demonstrate the consistency of this particle definition, by obtaining 
the same answers to this question as in [2, 6]. (Although the difference 
between these two questions is purely formal in this case, it becomes vital 
when considering gravitational fields, which cannot simply be switched 
off whilst making measurements.) The finite time behaviour of particle 
creation in the 'adiabatic' background is also studied here, providing a 
consistent account of the time-dependence of this creation process. We 
also consider the time-development of initial states containing particles. 

Section 2 is a short review of the approach to pair creation described in 
[13, 14]. It includes the particle definition appropriate for inertial observers, 
and a brief derivation of the general S-Matrix element and expectation value 
in terms of 'Bogoliubov coefficients'. Section 3 is devoted to calculating the 
finite time Bogoliubov coefficients in an spatially uniform, time-varying 
electric field in terms of the solution of a simple linear, second order equa- 
tion. In Section 4 we solve this equation in the case of a constant electric 
field, and use the result to find the number density of created pairs, and 
also their current density, after a finite evolution time. Section 5 presents a 
similar treatment for the 'adiabatic' electric field E{t) = (0,0,£^sech^(^)). 
Conclusions are set out in Section 6. 

2. BACKGROUND 

Let H = L^(IR^)^ (the space of finite- norm spinor- valued fimctions of 
space) be the space of 'first quantized' states. The antisymmetric Fock 
Hilbert space over H, denoted J^/^CH) = (B^^qA^-H, is the space of super- 
positions of all possible Slater determinants of first quantized states. The 
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time-dependent 'first quantized' Hamiltonian Hi{t) iH^H is defined by: 

Hi{t)^P{x,t)=^-ijykVk+mj°^i;ix,t) (1) 

where ak = 7^7^^ = 7fe7o and Vfj,ip{x) = dfj,ip{x) +ieAfj,{x)ip{x). The Dirac 
equation can be written in terms of this Hamiltonian as i'VQijj[x,t) = 
Hi{t)'ip{x,t). Its expectation value is the spatial integral of the '00-component' 
of the energy-momentum tensor: 

Jt=to 

We can use Hi{to) to split Ti into H = Ti~^{to)(B'H~ (to) where: 

'H~^{to) is the span of the positive spectrum of -ffi(to) 
'H~{to) is the span of the negative spectrum of -ffi(to) 

W+(to) is the set of all positive energy states, and H~{to) is the set of all 
negative energy states, as defined in the background ^^(fo) at time to. In 
the Canonical formalism this definition would correspond to 'Hamiltonian 
diagonalisation' [15, 16, 17, 18] of the second-quantized Hamiltonian 

H{to)= I d'xT°o{kx)), 

where i^i^) is the field operator. In the past, Hamiltonian diagonalisation 
has been criticised [19] for its reliance on an apparently arbitrary choice of 
hypcrsurfacc T,t„ , (given by t = to here) and time-translation vector field 
/c^ (= 6q here) . This arbitrariness has been a particularly difhcult problem 
when considering QFT in a gravitational background, but the difficulty can 
be satisfactorily resolved [13, 14] by specifying Yi-m a-nd in terms of the 
worldline of the observer (or particle detector). The choice t = to, k^ = 6Q 
then emerges uniquely as appropriate to an inertial observer stationary in 
the frame of the electric field. When restricted to this frame, this definition 
is also similar to that of Fradkin et al. [21] or Greiner et al. [22]. However, in 
[21, 22] the Hamiltonian is taken to be H' {t) = H\{t)+eA^i^{x,t) so that the 
Dirac equation becomes idiip{x,t) = H' (t)')p{x,t). Sriramkumar and Pad- 
manabhan [12] pointed out that this leads to a particle definition which is 
gauge dependent. This problem, which has plagued previous 'Bogoliubov 
coefficient' approaches to particle creation, is usually patched up by appeal- 
ing to the tunnelling interpretation: if a gauge A = {A{x), 0,0,0) is chosen. 
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particle creation is described in terms of particles 'tunnelling through the 
barrier that separates particle/antiparticle states' (see [3, 22] or [5] for a 
description of these methods). This leads to two theories, each succeeding 
where the other fails [12]. This is not the case here, however, since we use 
eigenstates of the gauge-covariant Hamiltonian Hi{to) to define particle 
states. It will emerge that particle creation in the background re- 
quires either: 1) that A^{x)^0, so that A°{x) can act as a potential, w.r.t. 
which tunnelling may occur, or 2) A{x) = {A^{x),A'^{x),A^{x)) is time- 
dependent, so that eigenstates of Hi{t) mix during evolution. Both effects 
contribute to the calculation of Bogoliubov coefficients [13, 14], and ensure 
their gauge invariance. An immediate consequence is that if a gauge can be 
chosen such that A° = and A is time-independent, then particle creation 
will not occur. That is, time-independent magnetic fields cannot create 
particles. This result agrees with Schwinger's calculation [1] in terms of an 
effective Lagrangian, and siicc;cssfully resolves the inconsistency problems 
raised by Sriramkumar and Padmanabhan [12]. 

Having defined (to), we can now define the vacuum at time to, jvact^), 
to be the Slater determinant of any basis of TL^ {t^) (normalised so that 
(vaCfojvacto) = 1). This specifies jvacto) up to an arbitrary phase factor. It 
is the state in which all negative energy degrees of freedom are full, and 
hence is a concrete manifestation of the Dirac Sea. 

To illustrate this, let {wj^toiie/}, {fi,to;i€-f} be orthonormal bases for 
W+(to) and 'H~{to) respectively, where I is some index set, assumed count- 
able for convenience (the uncountable case introduces no complications). 
The vacuum at time to can be written as: 

|vacto)=wi,toAt;2,toA... (2) 

and is independent of the choice of basis for (to) (up to a phase factor) 
because of the complete antisymmetry of the Slater determinant. The 
vacuum at some time ti > to is: 

|vacti)=Ui,tiAw2,tiA... (3) 

By considering fermionic QFT 'in a background' we are ignoring direct 

electron-electron interactions, so that the evolved state \va.Cto{ti)) obtained 
by evolving jvacto) from time to to time ti is simply: 

I vacto (i 1 )) = -fi ,to (t 1 ) A V2,to (i 1 ) A . . . (4) 

where w;./„(ti) denotes the state obtained from Vi^to by evolution to time 
ti, and will not, in general, be contained in 'H~{t\). We will often refer 
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to this state as the 'evolved vacuum', although it is not in fact a vacuum 
state. 

Since the inner product of Slater determinants is simply the determi- 
nant of the 'first quantized' inner products, the vacuum-vacuum S-matrix 
element is simply: 



(vacti |vacto {h)) = det[{vi^ti \vj,to {h))] 



(5) 



vac— ^ vac ' 



The probability that |vact(, {ti)) will be vacuum at time ti is then 
|(vac(i |vacto(ti))p. Although once a particle interpretation is specified this 
result can be derived by a number of methods [21, 7, 23], wc believe that 
the derivation presented here (and in [13, 14]) is the shortest and clearest. 
Similarly, the general S-matrix element can be written as: 



(( 



:(-)''-^'det 



7in 



\missing/ 



R., ... ( ii '-i" ) 



\missing/ - 



en 



_\missing' 



( 31 ---U ~ 
\missingj - 

,. /3l-3n V 

\missingJ 



(6) 



{if m — n = m' — n' , and zero otherwise) where 



aijiti,to) = (Mi,ti|uj,to (ti)) 



ItjitiM) = {vt,tAuj,to{ti)) (7) 
eij{h,to) = {vi^ti\vj,to{ti)) (8) 



are the time- dependent Bogoliubov coefficients, {—Y is an unimportant 
sign convention, and K^^j^ j"')^ ) denotes a state comprised of m particles 
(in states Ui^^to ■ ■■'^im,to with ii <i2 < ■ ■ - im by convention), and n antipar- 
ticlcs (corresponding to the absence of states Vji,to---'Vjn,to)^ prepared at 
time tQ. 



Expectation Values 
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Given an operator Ai{t) : W — > W, we can define its physical extension 
^physW:^A(H)^^A(H) by: 

lphys(0 =^ifW - (vact|lij(t)|vact)T (9) 

where Ah : V'l A V'2 A • • • A V'jv ^ ^ V'l A • • • (Ai'^,) A V'i+i • • • A V'jv (10) 



The expectation value of Aphys{t) in the physical vacuum at time t, |vact) 
is zero by construction. Its expectation value in the 'evolved vacuum' 
Ivacjo {t)) is in general non-zero, and takes the form 

( vacto (t 1 ) I Aphys (i 1 ) I vacto (t i ) ) (11) 

N N 

= E(^^^,to(ii)l^i(ii)l^*,to(ii))-EK*il^i(ii)k^.ti) (12) 

i-l i=l 

= Trace(/3/3l' A++ -77t A~ +€/3tA+- +/3eU-+) (13) 
where we have defined: 

=i'^3M\Mil)\ukM) A'Jk={'U3M\Mil)\vkM) 
^jV = ("j,*il^i(*i)I^Mi> and A-+ jli(ti)|ufe,t J (14) 

= if A\ is Hermitian 

The operator A^phys(t) that represents the number of particles (includ- 
ing antiparticles) present at time t is the physical extension of Ni{t) = 
P^{t) — P~{t), where P'^it) ■.H^H'^{t) are the projection operators onto 
W±(t). It is easy to verify that ^phys(io)| (j;;^:.^)^^) = ('^+")l ( 
as required. The expectation value of iVphys(ii) in the 'evolved vacuum' 
|vact(,(ti)) is given (from (13)) by: 

N^acto itl) = (vaCto (t 1 ) I /^phys (t 1 ) | vact^ (t 1 ) ) 

= Trace(/3/3^-|-77^) (15) 

i 

where N^^^{ti) = {(3(3^ )u is the expectation value of the physical extension 

of I and represents the probability that the degree of freedom 

Ui^ti is occupied in jvactg (fi)), i.e. that particle i is present. -^j~tg(ii) = 
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(77^ )m is the expectation value of the physical extension — \vi^ti) {vi,ti \ ] it 
represents the probability that the degree of freedom Vi^ti is unoccupied in 
|vacto(ti)), i.e. that antiparticle i is present. From the unitarity of the 'first 

a{ti,to) P{ti,to) 
7(ti,to) e{ti,to) 

show that Trace(/3/3^) = Trace(77^), expressing charge conservation. 

The derivation of {Ft„{ti)\Aphys{ti)\Fto{ti)) for an arbitrary state |^to(ii)) 
is identical to the derivation of (13), and gives: 



quantized' evolution matrix Si{ti,to) ' 



, it is easy to 



m n 

= I]("H,to(*l)l^l(il)l«4,to(*l))-$Z(^^J.,to(*l)l^l(*l)l^^i.,to(*l)) 

k=l k=l 

+ {va.Cto{ti)\Aphys{ti)\va.Cto{ti)) (16) 
More details of this approach to fermionic QFT are given in [13, 14]. 



3. SPATIALLY UNIFORM ELECTRIC FIELDS 

A spatially uniform electric field in the z-direction can be represented 
by A^{t) ^ {0,0, 0,A{t)), so that E{t) ^ {0,0, A{t)) and B{t)=0. To find 
solutions to the Dirac equation in this background, we try separating ip{x) 
as 

^p{x,t)=iljp{t)e^P'' (17) 
where p-x = Y^l.^-^p''x''. Then V'p(*) must satisfy: 



'^^^S^ = m^oMt) +P±M*) + {p^+eA{t))a3ijp{t) (18) 

where pj^ =p^ai +p^CT2. By introducing a set of basis spinors 4'±±{p) 
satisfying: 

4'lx{p)(t^s'y{p) =Sss'S\x' for all p 

CT3'/>±a(p) =±(j)±\{p) for A = ± and all p 

locf>±x{p) = (f>T>. (P) for A = ± and all p 

P±(ps\{p) =sX\pj_\(p-s-\{p) for s = ±,A = ±, and allp 
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where p = {p^)^ + (p^)'^, we can expand tpp{t) as: 

Mt) = J2{^+x iP)h (t) + ct>-x{p)9x m (19) 

X 

The inner product of two solutions il)p{x,t),'tjjq{x,t) expanded as in (17) 
and (19) is given by: 

(4(a;,i)|V|(x,t)) = (27r)35(p-q)^(/l(t)/|(t)+ffi(t)ff^(t)) (20) 

A 

Substitution of (19) into (18) gives: 



i^-ip^ + eA{t))f±=mg±±\p^\g^ 

i^ + {p' + eA{t))g±=mf±T\pJh (21) 
Define now the dimensionless quantities 

T = mt Pz = — p=-^ and a(T) = — (22) 
mm m 

This means wc are measuring distances in multiples of the Compton wave- 
length ^ = — « 3.9 X 10"'^'^ metres, and times in multiples of « 

" 27r mc ' ^ 27rc 

1.3x 10~^^ seconds (where we have used the mass and charge of the elec- 
tron). Equations (21) become: 

i-^-{Pz + a{T))^f±{T)=g±{T)±pg^{T) 

'j^ + ^P^+''^^)))a±iT)=f±i^)^Ph(r) (23) 

from which we find that 



-^^yl+p' + {p.+a{T)Y+i—yx = Q (24) 



^ + (l+/ + (p. + a(r))^-zg)5. = 



(25) 



for both A. To construct a complete set of solutions to (18), let /(r) be 
any solution of (24). One solution of (18) can then be found by putting 
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/_(r) =0, /+(r) =/(r) and using (23) to get: 

1 p 

where /i(r) = {i-4- — + a(T)))/(r). Prom this we obtain 



(26) 



,ipx 



(27) 



where C= ((l+p^)|/p + |/ip)~2 is constant by virtue of conservation of 
the norm. Similarly, ■ip'p2{x,t) can be found by putting /+(t) =0 and 
/_(t) := /(t), and then using (23) to give g+{T) and g-ir) (in terms of 
i^p,3{Xjt) is found by putting (7_(t) = 0, g+{T) — J{t) (which solves 
equation (25)) and then using equations (23) to specify /±(t) in terms of 
/i(t), while -0^.4 (35, t) is similarly found by putting 9+ (r) — 0, ,g^(r) = /(r). 
This finally allows us to generate the complete set of solutions: 



V'p,i(a;,i) = r- M }'jjAx,t)-pi}'jjo{x,t)) 

Vi+p^ 

= C {p)-pcl>+- (p))/(t) + <^_+ {p)h{r)] e^P-^ 

V'p,2(a;,i) = —^^={p^pp^i{x,t) + 'il;p^2ix,t)) 
V^+P 

= C [(P0++ (p) + (p))/(r) + 0- (p)/i(r)] e^P^ 



V'p,3(a!,i) = C 
V'p,4(a;,i) = C 



t^±(^l±4±^Mr)+0-.(p)v/TT^/>) 

^^^^±±^fL|±^Mr) + 0--(p)v^/>) 



(28) 
jpx 

^ipx 



It is easy to verify that these solutions are mutually orthogonal and are 
each normahsed to (27r)^(5(0). The choice of ^p^i and tpp^2 here (in terms 
of tp'p,! and V'p,2) is for later convenience. 

3.1. Energy Eigenstates 

We now calculate H^ito) at each time to- We must first decompose the 

space of spinor-valucd functions ijj{x,tQ) (of space, at time to) in terms of 
the spectrum of the operator H{to) given by: 
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H{to) ■■ i>{x,to) -i'^akdktp{x,to) + eA{to)a3ilj{x,to) +m'yoi^{x,to) 



(29) 



Consider the action of H{to) on plane wave states of the form (17). It is 
easy to see that: 

H{to)'ipp{x,tn) = Ep^toi'pixM) if and only if 
mjoijp{to)+p_Ltpp{to) + {p^+eA{to))a3iJp{to)=Ep,t,ijp{to) (30) 

By substituting Vp.to =J2x{<P+>^iP)f^,to+'l^-\{p)9x,to} and using (22) 
we transform this into the eigenvalue problem: 



m 



Pz + a{T()) 1 

Pz+aiTo) -p 

1 -P -Pz-a{To) 
p 1 



P 
1 


-Pz-a{To) 





f+.to 




.f-,ta 




9+.to 




5-,io 



f+.to 
f-,to 

9+,to 
9-, to 



(31) 



This matrix squares to {m^ + \pj_\'^ + {p^ + eEA{to))'^)l4. Its eigenvalues 
are therefore Ep_i„ = -k^Jm? + p + [jfi ^ eE A(tQ)Y' = irnEr„ in terms 
of the dimensionless quantity Et = y/ 1 + {pz + a(T))^. We consider the 
following cases: 



Case 1: 



'P,to ■ 



-niEr 



We seek two independent solutions of (31) with -Ep,to ='mET^. So that 
these are similar to the ^p^i{x,t), we choose up^i^to{x) such that g+,to = 1 
and (?_^to = 0. In (31) this gives: 

I. r, (l}++{p)-P<i>+- jP) + {Ero -Pz- a(To))0-+(p) ijj.x f . 

Up,-i,to[x) = D — — (32) 

^To—Pz — a[To) 

where D = constant, and wc have used the identity {Et^ —pz —a(To))(Era + 
Pz+a{To)) — l+p^- By requiring Mp,i,to(a^) to have norm {2n)^6{0) we 
obtain!)- /^-o-P-«(-o) 



12 



C. E. DOLBY AND S. F. GULL 



Similarly, up^2,to{^) be found by putting = 0, g-^to = 1) ^^'^ then 
normalising to {2it)^S{0). This gives: 

/X <I>++{P) -P(t>+~ jp) + {Ero -Pz- a(To))0-+ (p) ip.x 

Wp,i,to (a;) = , e ^ (33) 

V 2^ro (^'ro - - a(7"0 ) ) 
/ N + '?^+- (P) + (^ro -Pz- a{To))<j)—{p) jjjx 

^2£;ro(^T-o-Pz -0(7-0)) 

Case 2: Ep^t^=-'mEr^ 
A similar procedure gives: 

\P) ip-X /orX 

V2£;ro(^ro+Pz + a(To)) 

- - ' ^+- (P) + (^^0 +1>^ +ai rom-{p) ^ip.x (36) 

^ V2i?.o(^^ro+Pz + a(ro)) 

These states are eigenstates of the gauge invariant 3-momentum operator 
—idk + eAk of momentum p/ = {p^,p'^,p^ + eA{to))- Hence the 1-electron 
state |(^''^)(^) = Mp,A,to ^ |vacto) physical (gauge invariant) momentum 
p/. The l-positron state \{p ^ ) = ivp a t(, ^ I'^^^to) (where i^p ^ removes 
the vp^x^to degree of freedom from |vaCfo) [13, 14]) has physical momentum 
-p. 

Suppose now that A{to) =Foo as to ±00. This holds for the constant 
electric field A{t) = Et {E < 0) and many other electric fields. In this case 
we have: 



X 



(^3Up,X-ooix) = Up^X-ooix) Cr3Vp^X-oo{x) = -Vp^X-oo{' 

CT3Wp,A,<x> (a;) = - wp,A,cx> (a;) ^3Vp,\,oo {x) = vp^x,oo {x) 

It follows that any electron (of charge e < 0) in this electric field, that has a 
finite momentum at finite time, will at late times be moving in the negative 
x'^ direction with infinite momentum (i.e. at speeds approaching that of 
light), whereas a positron with charge — e will be moving in the positive 
direction with infinite momentum. At very early times this situation is 
reversed. Clearly it is important to keep track of what constitutes parti- 
cles/antiparticles at different times, since not only does this change with 
time in an electric field, but here the spaces have completely reversed 
between early and late times; W^(— 00) =W^(oo)! This fact is also men- 
tioned in [7]. 
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3.2. Finite Time Transition Probabilities 

We now have the positive/negative energy eigenstates {\up,x-to) ,\vp^\-to)} 
at each time to, and the evolution operator: 



Ui{ti,to)-. 



(2^ 



^IV,»(ii)>(^p/,i(to)| 



1=1 



Hence we can write |wp,A;to(^i)) =C^i(ii>io)|Mp,A;to)> 
\vp,\-to{ti)) = U\{t\,tQ)\vp^\-to) and the time-dependent Bogoliubov coeffi- 
cients become 



Sp,X;q,cr{tl,to) 



o:p,x;q,a{ti,to) 0p,X;q,aiti,U)) 

lp,\;q,cr{tl,to) ep^X;q,aitl,to) 

{up,x,tAx)\Ui{ti,to)\uq.,j;to{x)) {up,x:tAx)\Ui{ti,to)\vq^^.to{x)} 

{'Vp,\;ti {x)\Ui{ti , to)\uq^a;to {x)) {vp,X;ti {x)\Ui{ti,to)\vq^a;to {x)) 

A 



Afp, A;p/,i (ti )Mq^a;p/,i (tq ) 



(37) 



where: 



Mp^X;P>,ii'^)- 



{up^X;t{x)\lpp,,i{x,t)) 
{Vp,X;t{x)\'>Pp'Ax,t)) 



(38) 



which is unitary for all time, since it describes the overlap between two 
orthonormal sets of states. 

We can now use (33) - (36) and (28) to calculate M(r) as 



Mp,A;p/,i(T) = {2nf6{p-p/) 



Bp,p,{T)l2 Ap^p^{T)l2 



(39) 
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where 



C 



^/2Er{Er-Pz-a{T)) 



{{E^-p^-a{T))h{T) + {lW)f{r)) 



(40) 
(41) 



Bp,p^ (r) ^ 



C 



^2Er{Er+p, + a{T)) 



{{Er+p, + a{T))h{T)-{l+p')f{T)) 



(42) 
(43) 



and I2 is the 2x2 identity matrix. From (39) we have 

5'p,A;qf,a(Tl,To) = {2TTf5{p-q) 

where: 



- 0p,Pz (n , To ) /2 ap,p^ (n , To ) h 



(44) 



Q!p,p, (n,To) = ^p,p^ (n)Ap,p, (to) + -Bp,p^ (Ti)-Bp,p^ (ro) (45) 
(n,To) = ^p,p^ {Ti)Bp^p^ (to) - Bp,p, (r^^p.p^ (to) (46) 

We can now use the formulae in Section 2.3 to extract from (44) what- 
ever information we wish about this system of electrons evolving in a back- 
ground E= (0,0, ^(t)), expressed in terms of an arbitrary solution to equa- 
tion (24). For example, suppose we start the system at some initial time to 
in the vacuum state Ivact^), and let this state evolve until some later time 
ti. Then, from (15), the expected number of particles per unit volume in 
this evolved vacuum (not including antiparticles) with spin A and physi- 
cal (gauge-invariant) 3-momentum p/= {p^ + eA{ti)) in a range (Pp 
about p/ is given by: 



p,\,to 



{tl)d^p 



V 



■■\l3p,p,{n,To)\^d'^p 



(47) 



(where V = {27r)^S{0)). The the same number of antiparticles is also created 
with physical momenta — p/. The total number of particles per unit volume 
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in this evolved vacuum (including antiparticles) is therefore 

(48) 

-.3 r°° 



= — 2" / dp^ / dpp |/3p,p,(ri,To)|^ 

The probability that this state is vacuum at time ti is determined from 
(5), using det(e^) = e''^*''''("*\ to give 

^^-2-/ dp^ / dpj3 log(l-|/3p,p^(Ti,ro)P; 

(49) 

From (33) - (36), the total current density in this evolved vacuum is (13) 
and Ji,k = eak) 

' dp^ -F^{{Pz + a{Ti))\fip.^pATi,To)\ 



- ^1+^ 5i(/3p,p,(n,ro)ap,j,,(Ti,To))} (50) 

where Jvac,i,to(*i) = = •^vac,2,to(*i) by symmetry. This expression will 
generally contain a logarithmic divergence proportional to a(T), which must 
be isolated and absorbed into a renormalisation of the coupling constant e. 
In Section 4, a(T) =0 so that no renormalisation is necessary. For most of 
Section 5 we use a cut-off in the integrals such that the cut-off-dependent 
renormalisation term is small, and can be ignored. 

Non- Vacuum Initial Conditions 

Consider an initial state at time to consisting of N particles having mo- 
menta Pi,...ppf and spins Ai,...AAr. From equations (16) and (44), the 
expected number of out particles at some later time ti is: 

AT 

AfouMo(ii)=I](|ap.,P...(n,ro)|'-|/3p„p,.,(Ti,To)n + 7Vvac,to(ii) (51) 



which is less than A'vac,to(*i) + -^ due to Pauli exclusion (recall that 

I "p,P. ("^1 > To ) P + 1 /3p,p, (n , To ) I ^ = 1 ) . 
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The probability that the out state contains exactly N particles is: 

?^iv^^.ofti)= ^/'r"°^^'7^„ (52) 
ni=iFPi,Pz,*(n,To)|^ 

This expression highlights the well-known fact that particle creation is less 
likely in a state that already contains particles. The current density of this 
evolved state is: 

ih) (53) 

i=l 

where from (44) 

Up,\,to ix,ti) = Ctp^p. (n , To)up, A,ti (x) - Pp^p^ (ti , To A,ti (x) 

SO that 

J{up,x,to{x,ti)) = -^{{l-2\Pp,p^{TuTo)\^){p, + a{n)) 

+ 2 ^l+p^^{ap,p^ in, To )f3p,p^ (ri , To ) ) } (54) 
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Simplification 

Although M{t) depends on the initial conditions chosen for /(r) in (24), 
the elements of S{ti,to) do not depend on them. We can, however, by an 
appropriate choice of initial conditions, simplify the expression for S(ti,tQ). 
Impose temporarily the conditions /(tq) = 1 and /i(to) = Etq —Pz — a(Tb)- 
(This second condition is equivalent i^\r=To =^'to/(t'o)i revealing a con- 
nection between our particle definition and that motivated by 'instanta- 
neous frequency'.) These conditions yield Ap^p^ (tq) = 1 and Bp p^ (tq) = 0. 
Hence we can identify ipp,i{t) = up^i^toit) and ■!/'i?,i+2 (0 = i'p,i,fo(i) for 
i = l,2, so that Q!p,p;,(Ti,To) = Ap,p^(ri) and /3p,p^{Ti,To) = -Bp^p^{Ti) in 
this case. The current density in the evolved vacuum then simplifies to: 



Jvac,3,to{ti) = / dpz dpp i 

Ero+Pz + a{To)) 



-|/(n)| 



Er^+Pz+a{Tl)) 

(55) 



Eti 



and J{up.^\.^to{Xtti)) becomes: 

J{up„Xutoix,ti)) = e<^ — |/(ri)| -1 

These simplifications will be used in Section 5. 



(56) 



4. CONSTANT ELECTRIC FIELD 

For constant electric field we have o(r) = ^t= aor. Since e is negative, 
positive 0,0 corresponds to an electric field in the — z direction. By writing 
ao = — we identify the 'natural' electric field strength Ec = « 1.3 x 
10^® V m~^ (where we have substituted the mass and absolute charge 
of the electron). This corresponds to an electromagnetic energy density of 
approximately 137 electron rest mass energies per Compton volume {-^) , 
and is too large to be produced on macroscopic scales in the laboratory, 
although astrophysical examples have been considered [24, 25]. Also, fields 
of this strength would be present on a microscopic scale in superheavy 
nuclei with charges greater than 137|e|. 

Equation (24) can now be written as: 

-^ + {l+p'^ + {pz + aoTf+iao)f = (57) 
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Comparison with equation 8.2(1) of [26] reveals that solutions of this equa- 
tion are linear combinations of Dv^-i{±^^^{pz+aQT)) and D-^^ (±-ii= (pj,+ 

aor)), where Vp = 2oo ^ ^^'^ D^{z) are parabolic cylinder functions. We 
choose f{T) = D^^-x{^^^{pz + aQT)), which gives (after exploiting proper- 
ties of the parabolic cylinder functions) 



h{T) = {l+i)D,^ 



1 \ e 

{Pz+aoT) I and C - 



8aQ 



^/2 



The first point about these expressions is that /(r) (and hence App^{T) 
and -Bp,p^ (r)) depends only on pz and r through the combination Pz + clqt, 
which is the physical ^-momentum (in units of mass) at time r of the state 
labelled by the subscript p. This allows us to write A.p,.p^{Ti) = Appout(O) 
where p°"*=p2 + aoTi, so that the evolution matrix can be written as: 



ap_j,out(0,-T)/2 /3p,p;"t(0,-T)/2 

-/3p_pout(0,-T)/2 Q:p,pout(0,-r)/2 



and depends only on the evolution time t = ti — tq. This follows from the 
time-translation invariance of the constant electric field. 

We now recognise that the physically relevant infinite-time scattering 
matrix is: 



lim Sp,x;qA^-T) = {2^)H{p-q) 



^(0)/2 



j4„ 



Bp^pout (0)/2 Ap^pout (0)/2 



-_Sp,p-t(0)/2 



(58) 



where we have used asymptotic properties of pout (t) and Bp pout (t) and 
discarded unimportant phase factors. The expected number of particles 
per unit volume in the infinitely evolved vacuum is now 



Np,x{<x) 
V 



= \B., 



p,pOut (0)1 
-7r(l+p^) 



(59) 



4ai 



AEoiEo+pr') 

-{l+p^)D,,, 



{Eo+p°^'){l + i)D,^,+,2, 

2-1 „„* 



Pz 



(60) 



Curves of this expression for various values of hq arc shown in Figure 1. 
We see that Np^xi'^) is effectively zero for negative and approaches 

e "0 for large . in agreement with [11]. The result is plausible since 
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FIG. 1(A). Number density in an in- 
finitely evolved vacuum as a function of 

p°"*, for ao = 1 and p = (top curve), 0.5 
(middle curve), and 1 (bottom curve). 




FIG. 1(B). Number density in an in- 
finitely evolved vacuum as a function of p 

for ao = 1 and p°"* = — 1 (bottom curve) , 
(outermost curve), 1 (top curve) and 5 

(innermost curve). A plot of e "o is 
also included, which coincides almost ex- 
actly with the curve for p°"* = 5. 



any particle, once created, will be accelerated to ever increasing by the 
electric field. The number density of created antiparticles is identical, but 
with the sign of reversed. 

A common misconception is that *^'y = e "o independent or 

[2] . This misconception stems from taking ti ^ oo and tq — > — oo when 
defining 'in' and 'out' modes, inadvertently taking the limit as oo. 
A more relevant result, which allows Schwinger's formula to be derived in 
the large time limit, is described below. 

Figure 2 shows the expected number density A^p.a(t) in the evolved 
vacuum for ao = 1 and p—l, as a function of after an evolution time 
T = 7 (right), 4, 2 and 0.5 (single peak). A graph for r = oo is also included, 
which is seen to fit the left peak of the t = 7 line, and becomes constant 
at large p™*. The finite time behaviour demonstrated here agrees with 
that presented in [6]. Figure 2 has an interesting physical interpretation. 
At early times, particles are created mainly with close to zero (Figure 
1(B) shows that small p is also favoured), as expected, since these pairs 
require the least creation energy. At later times, the right peak represents 
the same early-time particles (which have been accelerated by the electric 
field), and the left peak represents recently created pairs, again at low 
momenta. Figure 3(A) shows '^y^P for 00 = 5, 1 and 2. 
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FIG. 2. Number density in the evolved vacuum, as a function of pj"*) ^""^ oo = li P = Ij 

and T = 7 (right), 4, 2 and 0.5 (single peak). A graph for t = oo is included, which is 
seen to fit the left peak of the t = 7 line and becomes constant at large . 
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FIG. 3(A). Total number of particles in a 

Compton volume of the evolved vac- 

uum as a function of t for ao = ^ (bottom 
curve), 1 and 2 (top curve). 




FIG. 3(B). Vacuum- vacuum transition 
probability for ao = 1 and V = 4QTz^m~^. 
The initial burst of pair creation reduces 
the exact transition probability (bottom 
curve) below Schwinger's prediction (top 
curve) at early times. 



From Figure 2 we see that for. large r, Np^x{T) can be well approximated 



by 



(61) 
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where 0{p°^^) is the Heavyside step function. This allows us to approxi- 
mate: 

j Np^x{T)dp°^'^aoTe=^^^ (62) 
from which we obtain the total number of pairs created: 

Ntot (r) w — 4— e for large t=— (63) 

This analysis correctly predicts the final gradient of the graphs in Figure 
3(A), and is the result most commonly referred to in the literature [6, 
2, 22]. However, we see in Figure 3(A) that for \E\ <^Ec the dominant 
contribution to the total number of produced pairs comes not from this 
final gradient, but from the initial burst of pair creation. Even for \E\ = 
\Ec, approximately 0.44 particles arc created in a volume in the first 
10~^° seconds, with only ~ 0.005 particles every 10~^° seconds thereafter. 
For \E\ = approximately 0.014 particles are created in a volume in 
the first 2.5 X 10^^" seconds, with only 8.6 x 10~^^ particles created every 
2.5 X 10^^° seconds thereafter. This difference becomes exponentially more 
pronounced as the ratio |-^| increases. 

Similarly, substitution of (61) into (49) together with the relation 

log(l — e "0 ) 

, gives 

^vac— >vac(7') ~ Cxp 

in agreement with Schwinger's original calculation [1] (also in [27], [2], [6]). 
In Figure 3(B) we show the exact vacuum-vacuum transition probability, 
and the large time approximation of (64) for ao = 1 (we have taken V = 
AOtt^ for convenience). 

Vacuum Current, and Induced Electric Field 

Figure 4 shows the integrand in Jvac,3,fo (^i) from (50) as a function of 
for various evolution times, and includes a graph of the ]/3]^ contribution. 

Wc recognise that at late times the dominant contribution to the inte- 
gral comes from the ]/3]^ term; the other term generates rapid oscillations 




-{eE)^^ 1 



4^3 



E; 



for r large. 



(64) 
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-0.02 






(A). At = V^ 



(B). At = 7ViE 




(C). At = oo 

FIG. 4. Current density in an evolved vacuum after time At = V eE (A), 7\/eE (B) 
and oo (C), as a function of p^, for ao = 1 and p= 1. The contribution is also shown 
(darker curve). In (B) and (C) this term provides the dominant contribution, with the 
5)i(/3a) term producing rapid oscillations about it. 



about this contribution. The term is the 'classical' contribution of the 
form 'sum over created particles of the current of each particle', while the 
oscillatory 3?(/?q;) term represents quantum interference effects. These os- 
cillations are on a time scale tg„ ^ = ^ and a length scale of A^-^ Ac, 
as are the smaller, more strongly damped oscillations already present in 
the |/3p term. Oscillations on the time scale of tqy, can be identified with 
those mentioned in [11] and [28] as significant in quantum decoherence and 
effective dissipation processes. Figure 5 shows the total vacuum current 
density as a hmction of r, together with the |/3|^ contribution, and the 
approximation that results from ignoring the 3?(/3q;) terms and using (61). 

The [/3|^ contribution shows good agreement with Figure 5 of Kluger, 
Mottola et al. [11] (though they use a different scale). Wc also sec that, al- 
though the 'semiclassical' |/3]^ term dominates after the first 10"^" seconds, 
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FIG. 5. Total current in a Compton volume of an evolved vacuum as a function 

of T. The contribution is also shown, providing an accurate approximation except 
at very early times. The late time approximation of (65) is also included, which has 
correct final gradient but misses the initial 'burst' of creation, as with the total number 
of created particles. 



the 'quantum intcrfcrwicx'' term is significantly more important in this first 
instant. Finally, although we have Ji,vac(T) = = J2.vac(''') by symmetry, 
the integrands of Ji,vac(''') and J2,vac(''") still display rapid oscillatory de- 
pendence on 

The large-time approximation of J3,vac(''"), obtained by ignoring the 
3?(a/3) term and using (61), is: 



■^S.vacMs 



ejeEfAt 



(65) 



This would generate an induced electric field in the positive z-direction: 



Eind{T)'' 



47r3 



(66) 



We have not taken back-reaction into account, so these results can only be 

IL I I^c I 

considered realistic at times when Eind ^ |-E'|) implying At <C V^^^TjJ^- 
For small electric fields (where eE<^m?) this approximation is valid for 
an exponentially long time. It never holds for an infinite time, however, 
so calculations of infinite time pair creation effects without including back- 
reaction have little meaning. In fact, for \E\ Ec this approximation is 
valid only for At-^ 10~^^ seconds! 



Non- Vacuum initial conditions 

As an example of a non-vacnum initial condition, consider the state 
'Wp,_^,A,to(^) ^ |vacto(t)), which at time to contained one particle with mo- 
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mentum p^^. Recall from equation (54) that J{up^^,x^to{t)) does not rep- 
resent the current at time t of this particle. Rather, J{up^\^to{t)) + Jva,c{T) 
represents the expected current at time to + of this state, which can 
no longer be viewed as containing one particle. This current will contain 
contributions not only from the 1-particle component of the evolved state, 
but also from the entire evolved state. Notice that, from equation (50), 
Jvac{T) contains a factor V, while J{up .i.iaW) docs not. The vacuum 
contribution is that of a uniform current density pervading all space, and 
the contribution from J(up,A,to (*)) represents a single particle moving on 
this 'sea of current'. 




2468 10 2468 10 



(A) .\E\=Ec (B). |£| = ^ 

FIG. 6. J('ttp,A:to(3;,fo + ^ function of t, for p-^^^ = and \E\ = Ec (A) and 

(B). Also included is the classically expected current for a particle accelerated from 
rest in this electric field. As expected, the full result is closer to the classical result in 

(B) than in (A). 

Figure 6 shows the z-component of J{up^x^t„{x^t)) as a function of r, 
for Pijj = and \E\ =Ec (A), and j§ (B). Also included in these graphs 
is the classical current for a single particle accelerating from rest in the 
relevant electric field. For \E\ —E^ the classical result differs significantly 
from the quantum result, even when the vacuum contribution Jva.c{T) (from 
Figure 9) is included. This can be understood as follows. From equation 
(11) the vacuum current density is simply the sum over all modes of the 
difference between the two curves on these graphs (with the sign reversed, 
since it is negative energy modes that appear in (11)). We know from 
the curves in Figure 4 that this difference can be expressed as the sum of 
two contributions. The 3?(q!/3) term, which explains the 'zitterbewegung' 
oscillations in Figure 6, and the term, which explains the difference 
in height between the two graphs, and represents the probability that the 
particle falls into a hole in the Dirac Sea vacated by a newly freed vp degree 
of freedom. This interpretation is also supported by Figures (13), (14) and 
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(16) of the 'adiabatic' case. For \E\ = j^, for which the probabihty of 
particle creation is small, the quantum result agrees well with the classical 
prediction. The graph for \E\ = ^ is not included since the quantum result 
becomes indistinguishable from the classical prediction. Finally, if at time 
t we project out the 1-particle component of wp;_j,i,fo (0 ^ k^Ctg (i)) (for 
example if the wave-function 'collapsed' during a measurement of particle 
numbers) and measure the current of this projected state, we again obtain 
the classical 1-particle prediction. 



5. ADIABATIC ELECTRIC FIELD 

Consider now the spatially miiform electric field E{t) — Esech^{^). This 
is generated by: a(r) = aoptanh . Hence, (24) becomes 



+ 



ptanh [ — 
VP 



-iaosech 



/a = (67) 



To solve this equation we use the same procedure as in [2]. Define 



y^i(l + tanhQ) 



(68) 



so that we can rewrite (67) as: 



■{2 + ipAE)y}'^-u_{l + u+)g = Q (69) 



where l\E = E+^-E 

_oo is the change in energy caused by evolution from 
T = — oo to T = -|-oo in this electric field, Apz = 2aop is the corresponding 

change in z- momentum, and we have defined z/± = ^iAEiApz). It follows 
(see [26], pg 56) that a possible solution is the hypergeometric function 
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g{y) = F^iy-jl + v^-jl — ipE-^iy), from which we find that 
/(r) = F{u_ + 2/(r))e-'^— (1 + ) =^ 

Hr) ={^{l-y{T))F{l + i^_,l + iy+;l-ipE_^;y{T)) 

+ {E+^-p,-aQp)F{u_,l + u+;l-ipE_^;y{T))}e-''^—'^{l + e^)^^ 
C-2 =2E-^{E-^-p, + aop) 

(70) 

We have again made this choice of /(r) and /i(t) so that -Bp,p^ (r) ^ as 
T — » — oo and we can write S(ti, — oo) = Af (n) (dropping irrelevant phases) 
and use the simpHfications described at the end of Section 3 in this Hmit. 

Number density in an evolved vacuum 

Consider a state which at asymptotically early times contained no par- 
ticles. The expected number of particles (not including antiparticles) with 
physical momentum =p/= {p^ ,p^+eA{t)) in this state at some later 
time t is given by 

y -\^P,pAT^)\ 

As r ^ OO the asymptotic properties of hypergeometric functions [2, 29] 
imply that: 

Np^X-oo{oo) cosh(7rpAi?) — cosh(7rpAp2) 



V cosh(7rpA£') — cosh(7rp(i^+(x)+-B_oo)) 

Figure 7 shows ^ function of for ao = 1, p= 10 and 

IP I 

T = —5,0,5 and oo. The three curves in each figure are for ^^-^ =0,0.5 and 
1. These figures show how particle creation proceeds. Particles are created 
with small p^, and are then accelerated by the electric field to take their 
place in the final 'bulge'. The RHS of this final bulge comes from particles 
created at early times, while the LHS comprises particles that were created 
at late times and have undergone less acceleration since their creation. The 
height of the peak at the origin in Figure 7 depends on the strength of the 
electric field at that time, since this peak represents newly created particles. 
The corresponding distribution of antiparticles is obtained by reversing the 
sign of . 
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(C). T = 5 '"' (D). T = oo''" 

FIG 7. The particle density P->-.-^^ ^ ^jjg evolved 'in' vacuum as a function of 
for ao = ;^ = l, p=10, ■!^=0,0.5 and 1, and for T = mt = -5 (A), (B), 5 (C) 
and GO (D). 



Figure 8 shows the total number density of particles in a Compton volume 
(m~^) of the evolved vacuum, ^^^.--^t^) ^ function of t. This has been 
evaluated by Monte Carlo sampling over \px \ < 2, < 2 and over Ipzl < 8 
(which corresponds to — 3 < < 13) . This range is large enough to include 
the dominant contribution to the integral. The two curves are the result 
of independent sets of 100,000 samples, giving an idea of the numerical 
accuracy. 

Figure 9 is a repeat of Figure 7 but with oq = = 10 and p=l. Particle 
creation now occurs over a much shorter period of time, so the distinc- 
tion between the peak at the origin and the late-time bulge is much less 
pronounced. 

To understand this peak at the origin in more detail, examine Figure 

10. Here we show ^ as a function of r, for ao = ^ = 1, ^-^^ = 0, 

p = 5, and forp^ = -6 (A), -3 (B), (C), 3 (D), 6 (E) and ±40 (F^From 

(B), (C) and (D), the most significant change in the occupation number of 
any mode occurs at the time when the covariant momentum of that mode 
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FIG. 8. Monte-Carlo calculation of the total number density of particles in a Comp- 
ton volume for ao = 1 and p = 10. 



passes through zero. That is, a mode is most hkely to be excited from 
the Dirac Sea when its momentum passes through zero in the frame of the 
electric field. If the mode does not pass through zero (as in (A) (E) and 

(F)) then, although can become non-zero during periods when 

a(T) changes so that E{t) is significant, this effect is small and transient. 
This transient effect is also influenced slightly by the momentum of the 
relevant mode, with the position of the peak shifted towards times when 
the momentum was smallest. However, when the relative change in the 

momentum is small, the peak is displaced much less, and more 
closely resembles E{t). 

It is useful to compare Figures such as 7(D) and 9(D) with Figure 11(A). 
In Figure 11(A) we see that when oq is large the bulge becomes much more 
steep-edged. For ao = 200 we see that (at least for \pj_ \ = 0) any mode for 
which passes through zero will almost certainly be created. Figure 
11(B) shows how this result changes upon increasing Modes with 

higher \p±\ require more energy for their creation, hence an increase in 
IPj^I decreases the number and smoothes the distribution of the created 
particles. 



PARTICLE CREATION IN ELECTRIC FIELDS 



29 




(C). r = .5 '"' (D).T = oo''" 

FIG. 9. The particle density ^'^y ^ 'n the evolved 'in' vacuum shown as a function 
of for ao = ^=W, p=l, ■!^ = 0,1 and 2, and for T = mt = -.5 (A), (B), .5 

(C) and GO (D). 



Figure 12 shows the particle density — in an evolved state that 

was vacuum at time tq = mt = —5 as a function of for ao = ^ = 1, 

J^J- = 0, and for various p and r. The curves further confirm the present 
interpretation of the particle creation process. All have a peak at the right 
hand side (even for t = oo), which comes from the initial burst of particle 
creation, much as in the constant electric field case. The size of this burst is 
related to the strength of the electric field at time to, such that the graphs 
in 12(D) are the reflection of — about — a^p- Figures (A), (B) 

and (C) are qualitatively similar to the curves in Figure 2, especially for 
p = 20. This is expected, since a(T) does not differ significantly from ao for 
p = 20 -5<T<10. 

Current in an evolved vacuum 

Having seen how the process of particle; creation takes place, the pre- 
dicted current in the evolved vacuum provides no further surprises. 
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-0.002L 

(A), pz = —6. The covariant z- 




(B). Pz = —3. The covariant z- 
momentum of this mode is always nega- momcntuiu of this mode passes through 
tive. zero at t~3.47. 



0.4 




-10 -5 5 10 15 20 

-0.1 



(C). Pz = 0. The covariant z- 
momentum of tliis mode passes througli 
zero at t = 0. 

O.Olr 




-20 -15 -10 -5 5 10 

-0.002 

(E). Pz = 6. The covariant z- 
momentum of this mode is always positive. 




-10 -5 5 10 15 20 

-0.05 ■ 

(D). Pz = 3. The covariant z- 



momentum of this mode passes through 
zero at T Ri —3.47. 




-15 -10 -5 5 10 15 



(F). Pz = ±40. The position of 
the peak approaches zero as pz — » ±oo. 



FIG. 10. The particle density ^'^^ °° — in the evolved 'in' vacuum as a function of 
r, for 00 = ^ = 1, ■'^=0, p = 5, and forpz = -6 (A), -3 (B), (C), 3 (D), 6 (E) and 
±40 (F). 
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FIG. 11(A). The particle density 
— ^ in the infinitely evolved 'in' 

„3 

vacuum shown as a function of -aut for 



p = 10, 
and 200 



■ 0, and ao = -S- = 20, 100 



FIG. 11(B). The particle density 
— ^ infinitely evolved 'in' 

„3 

vacuum shown as a function of -ant for 
p=10, ao = ^ = 200, and -1^=0,5, and 
10. 





(B).r = 5 



(C). r = 10 (D). r = oo 

FIG. 12. The particle density P-y o '^^^ g^jj evolved state that was vacuum at time 
TO =mt = —5 shown as a function of , for an = -§- = 1, JSjJ- =0, p = 5,10 and 20, 
and for T = mt = (A), 5 (B), 10 (C) and oo (D). 



Figure 13 shows the integrand of "^'"'"'•'^y as a function of for 
aQ = ^ = 1, ^^J- =0, p= 5, and for various r; Figure 14 shows the same 
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(C). r = 5. (D).T = 20. 

FIG. 13. Current density in the evolved 'in' vacuum shown as a function of p°"* 
for ao = ^ = l, = 0, p = 5 and r = -5 (A), (B), 5 (C) and 20 (D). The 

contribution is also shown (the darker curve). The S)i(/3a) term at late times produces 
rapid oscillations about the \I3\'^ term, in contrast to early times (as shown in A). 



quantity as a function of r for various pz. The contribution has also 
been included. As in Figure 4, we see that at late times the 3?(/3a) term 
simply provides rapid oscillations about the |/3p contribution. Also, the 
contribution to from the = 6 mode is very small in Figure 14. 

This is because the covariant momentum of this mode never passes through 
zero. 

Figure 15 shows the total current Jvac,3,-oo(t) in a Compton volume 

of the evolved 'in' vacuum, as a function of r, for oo = 1 and p = 5. This 
current has been evaluated by integrating Jvac,3,-oo(i) over \px \ < 2, < 2 
and over |pz| < 8 (which corresponds to — 3 <p°"' < 13). This range is large 
enough to include the dominant contribution to the integral, but ensures 
that the cut-off-dependent logarithmic divergence term is small enough to 
be ignored. By analysing the high-momentum behaviour of the integrand 
we can confirm the existence of a logarithmic divergence proportional to 
a(r), as expected. The ciuves are not as smooth as those of Figure 8 
because the rapid oscillations have not been completely averaged out. 
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(C).pz=3. (D).p^ = 6. 

FIG. 14. Current density in the evolved 'in' vacuum shown as a function of t for 
ao = ^ = 1, = 0, p = 5 and = -3 (A), (B), 3 (C) and 20 (D). The 

contribution is also shown (darker curve). 



Non- Vacuum initial conditions 

Consider now a state which at asymptotically early times contained ex- 
actly 1 particle, with momentum p^^. As in Section 4, the current of 
this state at later times contains two contributions, one from Jvac,3,-oo(i) 
and one from J{up,^_x^^oo{x,t)). The ^-component of J{up,^^\_^oc{x,t)) 
is shown in Figure 16 as a function of r, for ao = 1, p = 5, \p™\ = 0, and 
p™ = —8,—5,—2 and 1. Also included on each graph is the classically ex- 
pected current for a particle accelerated from rest in this electric field. For 
(A), (B) and (C) these curves arc qualitatively similar to Figure 6, the only 
difference being that t now runs from — oo. In Figure 16(D), for which the 
physical momentum of the mode never passes through zero, the quantum 
result is indistinguishable from the classical result. 

6. CONCLUSIONS 

Wc have applied the particle definition proposed in [13, 14] to the study 
of particle creation in spatially uniform electric fields. By incorporating 
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FIG. 15. Monte-Carlo calculation of the total current in a Compton volume 
for ao = 1 and p = 10. The two curves are independent sets of 100,000 samples, giving 
an idea of the numerical accuracy. 



the 'Bogoliubov coefficient' and 'tunnelling' approaches into a single con- 
sistent, gauge- invariant definition, we have resolved several problems raised 
by Sriramkumar et al. [12, 8]. We have demonstrated the utility of a time- 
dependent particle interpretation by presenting the timc-dcvclopmcnt of 
the particle creation process, concentrating on the cases of a time-invariant 
electric field and an 'adiabatic' electric field. For a time-invariant electric 
field we found, as expected for a first-order linear evolution equation, that 
the particle content after a finite amount of time T is the same [2, 6] as after 
evolution of a free 'in' vacuum in the presence of a background that was 
switched on only for a time T. For the adiabatic case we have given a coher- 
ent account of the timc-dcvclopmcnt of the particle production process, in 
which the particles are created with small momentum (in the frame of the 
electric field), and are then accelerated by the electric field to fill out the 
'bulge' of created particles predicted by asymptotic calculations [2]. The 
current in the evolved vacuum is consistent with this picture, although the 
di{(3a) term corresponds to oscillations around the term and contains 
the well-known logarithmic divergence, proportional to a{t). We have also 
considered an initial state with one; particle, and described how this state 
evolves as the sum of two contributions: the 'sea of current' produced by 
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and p'/" = —8 (A), —5 (B), —2 (C) and 1 (D). Also included is the current expected 
classically for a particle accelerated from rest in this electric field. 



the evolved vacuum, and the extra current arising from the initial particle 
state. 

The repeated mention of 'small momentum' throughout this work may 
cause concern about the Lorentz covariance of our procedure. In a constant 
electric field the class of observers moving at constant velocity in the z- 
direction all measure the same electric field in their rest frame. Can they 
really all measure the particles to be created with low in their rest 
frames? Of course they do not all measure this. The reason is hidden in 
the initial conditions. Insistence that the state be vacuum at time t = tQ 
effectively chooses a fixed reference frame O. An observer in another frame 
O' would find no time t' for which the state prepared as above contains no 
particles. Consequently, no contradiction arises. 

We hope this work, along with the formalism presented in [13, 14], has 
shown the computational and conceptual value of working with a concrete 
representation of the Dirac Sea. We strongly support Jackiw's claim [30] 
that "physical consequences can be drawn from Dirac's construction" . 
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